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A numerical study of the quantum double pendulum is conducted. A suitable quantum scaling 

is found which allows to have as the only parameters the ratios of the lengths and masses of the 

two pendula and a (quantum) gravity parameter containing Planck's constant. Comparison with 

classical and semiclassical results is used to understand the behaviour of the energy curves of the 

- - . levels, to define regimes in terms of the gravity parameter, and to classify the (resonant) interactions 
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I. INTRODUCTION 

In 1992, J. Ford suggested the double pendulum as a suitable system on which to experimentally test the ability of 
quantum mechanics to describe classically chaotic systems [1]: as a simple spatially bounded autonomous Hamiltonian 
^**^ ■ system presenting -in its classical version- transition to chaos it appeared ideal to test the consequences of the lack of 
chaos in "eigenfunctions, eigenvalues and time evolution" of such quantum systems. Since then, our understanding of 
the signatures of classical chaos in quantum mechanics has greatly improved [2] ; in particular arguments have been 
H [ advanced to the effect that continuous observation of the quantum system would bring its behaviour back to the 
classical one [3] . Still the above mentioned characteristics of the double pendulum make it a candidate for a thorough 
quantum study: it is far simpler than other chaotic autonomous systems that have been studied in recent years, like 
quantum stadia [4], whith their relevance to quantum dot technology, or hydrogen in strong magnetic fields [5] and 
helium [6], which moreover are not spatially bounded; the quantum kicked rotor [7] and hydrogen in monochromatic 
^ , microwave fields [8] have been the focus of much study and have almost become a paradigm of "quantum chaos" , but 
they are non-autonomous. 

The relative freedom with which mass and lenght ratios can be changed in the double pendulum would moreover 
allow us to study how the classical transition to chaos is reflected in its quantum counterpart when varying these 
(^ ] parameters and this may be illustrative of various "quenchings" of chaos in systems like for example helium itself: 
^^ the absence of observed chaos in helium appears to be a consequence of the possibility of an adiabatic seperation 

^^ in hyperspherical coordinates [6] so that there are quite good quantum numbers right up to the present precision 
^H of observations. Different mass ratios would probably break this quasi-symmetry; but changing the electron masses 

• ^ would just result in making the system even more complicated and not accessible to direct experimental testing. 

r-| , On the other hand, even though the classical double pendulum has often been used as an example of autonomous 

chaotic system [9,10], doubts about the proper ordering of the non-commuting operators in its kinetic energy [11] 
have up to now made its quantum counterpart not palatable to extensive analysis. The aim of my paper is to propose 
a reasonable ordering giving the proper behaviours for zero gravity (the behaviour in the high gravity limit does not 
?H ' depend on the ordering), and to explore the properties of the eigenfunctions of the resulting Hamiltonian. To my 
knowledge this is the first time that the quantum dynamics of the double pendulum is analyzed at all values of gravity 
and not only in the high gravity limit where it reduces to the trivial case of two coupled harmonic oscillators. 

The present paper is thus organized: section II presents the system, both classical and quantum; section III in- 
troduces and discusses the quantum numerical methods used. Finally, section IV discusses the quantum system 
properties; first those that can be obtained by semiclassical methods, then those obtained from the quantum simu- 
lations, focusing on the different quantum behaviour in the three classical regimes of low gravity (classical regular 
motion in most of the phase space), medium gravity (classical "global chaos" regime) and high gravity (regular 
"coupled oscillators" regime). 

II. THE MODEL: CLASSICAL AND QUANTUM HAMILTONIANS 

An ideal double pendulum is shown in fig. 1; li and I2 are the lengths of the two pendula and Mi and M2 their 
masses. Introducing the ratios I = h/h and /i — M2/M1 and the momenta Li (total angular momentum of the 
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system) and L2 conjugated to the two angles (pi and (p2, the classical Hamiltonian reads [9]: 
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It is possible [9] to scale the system via the adoption of new adimensional variables: time r = ty^E/{Milf) and 
momenta A^ = Li/y/EMilf, where E is the constant total energy. The scaled Hamiltonian h = H/E then always 
equals 1 and depends only on the gravity parameter 7 — gMili/E and the ratios I and ji: 
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Classical results are presented through Poincare sufaces of section (SOS); in the present paper I shall only consider 
sufaces of section in the (p2 = 0, (^2 > 0, {Ai, (pi} plane of which several examples are given in Fig. 8 of Rcf. [9]. 

For 7 = the total angular momentum Ai is conserved, it is therefore one of the two actions of the system and 
the SOS consists of horizontal lines; motion in (pi is always a rotation, but -as ipi is not the angle associated to that 
action- tpi is not constant. The second action I2 has instead to be calculated numerically [9]. The uniformity of the 
SOS also hides the two different kinds of motion in (^2 = 0: rotation for |Ai| < ^2[1 + ^(1 — ly], and hbration for 

V2[i+Mi-o"] < iAii < .yw+J^o^^Wi- 

For 7 7^ we shall here only note two facts. One is the vertical asymmetry of the SOS which is due to the fact that 
-following Poincare prescription- only the orbits crossing it with ip2 > sue shown; the SOS for ip2 < is perfectly 
symmetric to it. The other is that for 7 > [2(1 + ^ + /^O]^^ rotation in (pi is no longer possible and pi is limited 
between ±arccos(l — [7(1 + n + fil)]~^). 

In quantum mechanics we cannot use the classical scaling; we instead multiply the Hamiltonian by 2Mil\/h so as 
to have as sole parameters the scaled (adimensional) gravity 7 — 2Mll\g/fi^ and again the two ratios / and /i; the 
adimensional scaled energy will be indicated as _E = 2Mil\E/h^ ^ the time as t — th/2Milf, and the adimensional 

scaled momentum operators as Li = Li/h where the quantum momentum operators are defined in the usual way: 
Li — —ihd/dipi. The symmetrized quantum Hamiltonian then reads 
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Infinitely many other symmctrizations of the last two kinetic terms are possible [12]; our choice has been dictated by 
the physical argument that for zero gravity the ground state is completely delocalizcd [13]; its energy must therefore 
be zero. As we shall see, the chosen symmetrization guarantees that this be the case, even for the truncated basis 
sets we have to use for our numerical simulations. 

Classical and quantum adimensional scaled parameters and variables are related thus: 
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Eq. (4) means that -for given values of I and /i- the energy levels corresponding to the same classical situation, as 
described by the SOS at a given classical gravity parameter 7, are to be found on the {7, E} plane along the straight 
line E — 7/7. 

The classical limit is obtained for 7 ^ 00 along such a line; this means having Mi, ^1 ^ 00 while keeping all three 
the classical parameters Z,/i, and 7 constant. The same result is obtained with the usual "unphysical" limit 7i — > 0. 
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III. NUMERICAL METHODS 

A. The Projection of the Double Pendulum on a Rotor Basis 

To numerically calculate the energy levels of the Haniiltonian (3) we project it on the basis $mim2 given by the 
tensor product of the bases for two free rotors: 
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and then we diagonalize the finite matrix obtained by truncating the basis at suitables values of the indices mi and 



TO2 [14]. The matrix elements (m[m2\H\mim2) = Jq /g rf<pid(^2<&J^',„/ -ff^^mima are 
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where we have separated the kinetic energy term K and the potential energy one U. 

It can immediately be seen that the kinetic energy matrix elements between the state mi — m2 — and every basis 
state (including itself) are zero, thus giving zero as an eigenvalue of the system for zero gravity. 

When running a simulation on a trucated basis, it is important to evaluate how many states are a good approxi- 
mation to those of the full problem. Our simulations for /j, ^ I ~ 1 use |mi|"*"^ = 25 (51 levels) and |r7Z2|""*^ — 18 
(37 levels) [15]; a test run at 7 = with a doubled basis set (81 x 49 is the best choice in this case) shows that 555 
levels (£;™°^ — 145) are practically identical to those calculated on the smaller basis, giving a total of about 30% of 
reliable states. With increasing 7 this number decreases, since the lack of interaction with the missing states at high 
energies will progressively make also the topmost reliable states unreliable. On the other hand, since the energy of all 
levels grows with 7, E"^"-^ grows at least as the slowest growing state, namely the ground state which grows as y^. 

We have performed other two reliability tests on (ordering independent) properties of the system: 

1) We confronted the Z = ^ = 1, 7 = level density with the theoretical value obtained from the third graph in Fig. 
3 following the procedure given in section IV A. As shown in Fig. 2, the agreement is very good for all the states we 
have found above to be reliable. 

2) Again for ji ~ I — 1, the dependence of the ground state energy from 7 approaches for high values of 7 the 

theoretical one from cq. (11): Eq^ = V 7(2 + a/2) ~ 1.848a/7. 



B. Husimi Functions 



For many years now, Husimi functions [16] have been widely used when comparing quantum and classical systems, 
as they allow to project quantum functions in phase space in a way that avoids the interpretation problems connected 
with Wigner functions. In Hilbert space, the coherent states to be used as coarse-graining functions are, for the 
cylindrical phase space of each of the two spatial variables i = 1,2, and apart from unnecessary constant phase terms 
[17], 
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Here the variables used are the adimensional quantum scaled ones, (pi and L, are the (quantum scaled) phase space 
coordinates of the center of the packet, and ai -the angular momentum width parameter- is a free parameter; good 
results are obtained when cr ~ 1 . The normalization is chosen so that the Husimi function of any single rotor eigenstate 
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is normalized to 1. 
The Husimi function for a double pendulum eigenstate |$) = 'S'mi.m2^rni.m2 |'™ii'7i2) will therefore be: 
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For comparison with the classical Poincare surfaces of section, we shall here calculate only the Husimi function (10) 
on the surface (Ai, (pi), where ip2 ^ Q and 
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IV. RESULTS 

Now that we have the necessary numerical tools, we can use them to explore the quantum behaviour in the three 
classical regimes which we encounter when increasing the classical gravity parameter 7 from to 00 [9]: regular 
motion in most of the phase space (7 ^ 0), "global chaos" regime, and regular "coupled oscillators" regime (7 ^ 
max(l/(2/x^), 1/[2(1 -I- /i)])). First though I shall make some general considerations which will help orientate us in the 
parameter space, and then pass to a detailed analysis of my numerical results. 



A. Level classification and densities 



Energy levels at a given 7 can be classified in three groups according to the character of the classical SOS they 
correspond to. Starting from the bottom we first have "coupled harmonic oscillators states" (high 7), then "chaotic" 
states (medium 7), and finally "free rotors" states (low 7). 

For low 7's only few states belong to the first two classes; for increasing 7, their number grows, but it remains finite 
for any finite value of 7; the number of free rotors states is instead infinite for every value of 7. 

For "coupled harmonic oscillators" states (in whose number is included the ground state), the energy levels -Eni,n2 
are given by the expression 
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and therefore grow as the square root of 7. In (11) the frequency factors are [9] 
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and the Maslov indices are both 1/2, as each of the two oscillators has two caustics (in this case the inversion points 
on the paths on which the actions are calculated) [16]. 

The number of levels under a given value E of the energy is therefore 



N'. 



so that, on one hand, the density of levels is: 



-2 
E /2 



~2 

E 



2701 q;2 27 



l + ^i 



dN 


E 


/1 + M 


1 


/l + M 


dE 


-7I 


/ I 


l\ 


/ / 



(12) 



(13) 



and, on the other, the number of levels under a given value 7 of the classical gravity parameter 7 grows linearly with 

7: 
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For "free rotors" states (7 = 0), it is more difficult to exactly evaluate the density of levels, as the second action 
cannot be calculated analytically. On the other hand some considerations can be made: while in "coupled harmonic 
oscillators" regime the energy is directly proportional to the scaled actions n^ (see eq. 11), in "free rotors" regime 
the energy is proportional to the actions squared (see eq. 5). Since the number of levels below a given energy is 
proportional to the product of the actions, in the latter case it is only linear in energy (as opposed to the quadratic 
dependence we have in the former case: see eq. 12) and the density of states is constant. This constant has to be 
evaluated numerically: to do it we start by plotting the classical scaled action I2 versus /i = Ai thus obtaining the 
constant energy curve at the classical scaled energy h (that by definition equals 1). We now note that increasing the 

energy the graph expands radially, the total number of states below a given energy E is therefore the number of 
couples of quantized actions that can fit in the area Ap swept by the constant energy curve in its growth from E — 

to E = E or -equivalently- the area swept by a radius connecting the points of the E = E curve to the origin. Care 
must be taken to calculate twice the areas swept both by the outer curves (corresponding to libration in the second 
angle ip2) and by the inner ones (rotation in (^2), as regions of phase space corresponding to different classes of motion 
have different quantum numbers [18]. We now recall the relationship eq. 5 between classical and quantum scaled 



actions that here reads rii = liy E/2, where rii = mi and ri2 sue the quantum numbers at energy E and li, i = 1, 2 

are again the classical scaled actions; the level density dN/dE therefore equals half of the area 2Ah as measured for 
the classical scaled energy h = 1. Examples are given in Fig. 3. 

In quantum scaled variables, the energy of the lowest states grows as ^/j; the energy of the highest states instead 
grows linearly with 7; the net result is a decrease in the density of states with 7 at any given energy E. On the other 
hand, when going to the classical limit (7 —^ 00, 7 = const), the density of states remains constant both because of 
the constant density at 7 = (for low 7) and because of eq. 13 (for high 7). 

In terms of physical variables the energy density instead grows going to the classical limit, due to the relationship 
^ = — j[^4^: but this contribution is a uniform scale one: it does not alter the level structure. 

B. Level Interaction and Relationship between Husimi Functions and Classical SOS for I = jj, = 1 

1, Level Structure and General Considerations 

Fig. 4 shows the energy curves for E up to 60 and 7 up to 10. As expected from our discussion in the previous 
section, we see that the lowest levels -which almost from the start are in the "coupled harmonic oscillators" regime- 
grow as -y/7; most of the other levels instead grow at first linearly with 7. This is a consequence of the adiabatic 
teorem for noninteracting levels: as the action is approximately constant, the growth in energy of these levels goes 
as the average potential energy < U >~ 7(1 + ^ + ^l) = 37. Thus the levels exibiting such behaviour must be (as 
confirmed by their Husimi functions) those associated at first with the surviving KAM tori [19] and then with the 
island chains with long recurrence times which are located at the highest values of |Ai| and at |Ai| <~ 0; in particular 
we shall see that the levels associated with the KAM tori at the highest values of |Ai| are very resistent and preserve 
the shape of their Husimi functions well into the "global chaos" region, where -if we look at any classical SOS- we see 
chaos almost everywhere. Only when the growth in energy of these levels slows down do the Husimi functions change. 

On the other hand, the levels associated with resonance islands grow from the start more slowly than the others 
[20]. Extremely noticeable are the groups of levels associated with the main resonance island just below the positive 
branch of the separatrix (Ai ~ \/2); those starting at ii^ = 1.5, 5.5, 12.0, 21.5, 33.8, and 48.2 are clearly visible in Fig. 
4. The 7 = energy of the lowest level of each group can be obtained only approximately from eq. (5): due to the 
energy and angular momentum discretization, we have that Ai in i? = 2 (mi/Ai) is not fixed, it varies -for scaled 
energies up to E = 140- between 0.81 and 1.21. We thus have for some mi values two groups of levels, both associated 
with the main resonance island: for mi = 2 these are at Ai = 1.20502 (E — 5.5) and 0.81379 (E ~ 12.0); for mi = 3, 
at Ai = 0.91544 {E = 21.5) and 0.73098 (E = 33.8); and for mi == 6, at Ai == 0.91653 {E ^ 85.7) and 0.81478 
{E ~ 108.4); note that, again from eq. (5), the lowest levels within each mi = const series are those at highest Ai. 

In Fig. 4, the two straight lines at 7 = 1/2/iZ = 1/2 (lower line) and at 7 = 0.11, where the last invariant torus 
disappears [9], (upper line) mark the region of classical "global chaos". This "global chaos" region appears darker 



than the rest of the {7, E} plane because it is there that the mi, —rrii degenerate levels significantly separate and the 
resulting high density of distint levels produces multilevel interactions of states corresponding to different classical 
resonances. These multilevel interactions correspond to the overlap of classical resonances and therefore are the 
quantum mechanical mark of "global chaos" [21]. 

Two-level interactions can be classified into two types: the first one is the splitting of degenerate or near-degenerate 
levels increasingly repelling each other with growing 7, which can be locally described by the Dcnikov model [22]: a 
two level Hamiltonian with constant diagonal terms and off-diagonal terms which depend on the perturbation (7 in the 
present case). The second type of level interaction is instead the avoided crossing, best described by the Landau-Zener 
model [23]: again a two level Hamiltonian, where now are the diagonal terms which depend on the perturbation while 
the off-diagonal terms are constant. 

Demkov-like level interactions are localized either at 7 = (interaction of near-degenerate resonance island levels) 
or in the "global chaos" region (breaking of the rrii, —mi degeneracy). Landau-Zener- like interactions are instead 
evident almost everywhere in Fig. 4. Still, the highest density of both splittings and avoided crossings is in the "global 
chaos" triangle, where the splitting of the toi, — toi degeneracy also induces a high number of avoided crossings. 

2. Husimi Functions at Zero Gravity 

At 7 = 0, the levels are all degenerate in pairs (except the rax = ones); the chosen basis set then decides how the 
probability is divided between the two states of each pair: looking at the two </?2 = 0, (Ai, ipi) planes (both (/32 > 
and ip2 < 0), the Husimi function of the mi > level of a degenerate pair has support on the upper half of the plane 
(Ai > 0) while the Husimi of the mi < one has support on the lower half. The size of the projection on each of 

the two planes instead depends on the underlying classical phase space structure: if \mi\ < y E we classically have 
rotation in ip2 ', one of the levels of the pair has therefore support either on the ip2 > plane or on the ip2 < one, 
with only a negligible tail on the other plane that becomes larger for states whose support is close to the sepatatrix; 

the other level vice versa. If instead |toi| > VE, classically we have libration in (^2: both the degenerate functions 
therefore have significant projections on both the planes. Levels with the same mi but with support one on ip2 > 0, 
the other on ip2 < 0, are not degenerate. 

3. Husimi functions at Low Classical Gravity Parameter 

For 7 '^> the pairs of levels are still essentially degenerate, but the interaction -though small- begins to mix states 
with different mi quantum numbers; in particular there is some flow of probability between states with opposite ttji's: 
again looking at two (quasi)-degenerate states, the tail of the mi > Husimi now also has a component in the lower 
half of the SOS; likewise, the tail of the mi < Husimi has a component in the upper half of the SOS. These tails 
-which for levels 410 and 411 are at 7 = 0.025 (7 = 2.3 • lO^**) already larger than the 7 = ones by five orders 
of magnitude- reflect the appearance of classical libration motion in ipi when 7 > 0; but, since they are due to a 
probability flow which happens via tunneling through the unbroken tori around Ai = 0, they remain small till the 
"global chaos" triangle is reached and those tori are broken. 

For 7 small enough that classical chaos is not yet global, and the resonance islands still cover most of the phase 
space (7 <~^ 0.1), the widest avoided crossings undergone by the main resonance states are with states having similar 
Ai but with the sign changed; no probability flow is visible around Ai = 0, again because of the unbroken KAM tori 
in that region; the flow takes instead place between the ip2 > husimi of one state and the ip2 < one of the other. 
Three examples arc given in Figs. 5, 6, and 7, together with the classical Poicare section for the parameters of the 
crossing; since the functions are symetric for (pi -^ —y^i, only half of the SOS is shown; the avoided crossings of the 
other states of each of the degenerate doublets are identical to those shown. 

The first example (Fig. 5) shows a very clean avoided crossing undergone by the third state of the group of states 
associated with the main classical resonance originating at i? ~ 86; the parameters of the crossing place it out of Fig. 
4, but it has been chosen because -being at higher energy- the Husimi functions are better localized in phase space 
and fiows at the crossing can be better recognized. Both levels being in the central region of the SOS (rotation in ip2 
for 7 = 0) the flow is completely perpendicular to the {Ai, ipi} plane: the two structures grow and fade but do not 
touch in the plane. 

The second example (Fig. 6) shows a similar avoided crossing for the first state of the same grouping; here part of 
the support of the second function (b) before the crossing is on the unstable fixed point of the principal resonance (it 
is therefore, at least in part a "scarred" state [24]); to the process already seen in the previous example is thus added 
a visible probability flow between the main resonance island and the scarred portion of the second state. 



Finally, the third example (Fig. 7) shows the avoided crossing of three levels, one of them being the first state of 
the group of states associated with the main classical resonance originating at E :^ 108.5. Of the other two states, 
one is its symmetric in Ai (apparently a scarred state centered on the unstable fixed point of the most noticeable 
resonance in the lower part of the SOS), the other one (the intervening state) is instead a mixture of a scarred state 
of the period two resonance above the principal one with an excited state of -again- the most noticeable resonance in 
the lower part of the SOS. Again, even though the support of a state can move to Ai > to Ai < when passing the 
avoided crossing, all visible flow in the {Ai, ifi} plane is among structures with the same sign of Ai. 

From our study of the Husimi functions accessible to my simulations in this regime, it appears that level interactions 
which can be described by the Demkov model (level splitting) take place in the {Ai, (pi} plane, while interactions to 
be described by the Landau-Zcncr one (avoided crossing) thake place perpendicular to it. The first part of the above 
statement could be expected from what wc have already seen -namely that the Demkov transitions are those mixing 
states with different values of mi- the second part is instead a consequence of the mixing of states with different 
values of n2 by the Landau-Zener transitions. 

4- Husimi Functions in the "Global Chaos" region 

As we have seen, the two processes responsible for the high density of avoided crossings in the "global chaos" triangle 
are the growth in energy of the low lying states and the splitting of the ±TOi degeneracy. Both these processes are 
not associated with the appearance of resonance islands in the {Ai, ipi} SOS as these latter are connected with the 

interaction of states with similar mi/vE ratios but different |?tii| and n2 values. 

At such 7 values the phase space is mostly taken by the chaotic sea and little remains of the classical resonant 
structures visible at lower 7 values but many quantum states still have Husimi functions peaked on their stable and 
unstable fixed points. On the other hand, due to the multiple level interactions we have already mentioned, only 
rarely the Husimi fuctions of states in this regime are peaked on single structures: even away from avoided crossings 
the support of most states covers several classical structures, resulting in rather complicated multi-peaked Husimi 
functions. A few examples are shown in Fig. 8 a through d. Fig. 8 e and / instead show another typical shape for 
Husimi functions in this regime: the probability is concentrated along the border of the accessible classical region. 
Here -with the exception of the main resonance island which disappears at 7 « 2.0- are the last island chains to be 
eaten up by the chaotic sea (at 7 ~ 1.2) and the first ones to appear (at 7 w 2.5) when with increasing 7 the phase 
space reverts to regular. 

5. Husimi Functions at High Classical Gravity Parameter: Coupled Oscillators Regime 

For arbitrarily small quantum gravity parameter the ground state is in the "global chaos" region and its Husimi 
function extends on the whole SOS; when gravity is increased it soon leaves the "global chaos" region (at 7 == 0.075, 
7 is already bigger than 1/3) and at first concentrates on the fixed point just below the center of the 7=1 SOS in 
Fig. 8 of Ref. [9]; but when -at 7 « 8- it enters the coupled oscillators regime (7 ~> 2) where the SOS consists of 
concentric curves, it splits in two peaks located at Ai = and ipi close to the extreme values ± arccos [1 — (87)"^] 
(inversion points of the classical orbit). This reflects what can be observed in Fig. 8 of Ref. [9]: when with increasing 
7 the system leaves the "global chaos" region, the flrst regular structures to appear are at the rim of the SOS. For 
7 ^> 2 both the extreme values of (pi and the positions in ipi of the peaks decrease as 1/^/7 (or -equivalently- 7^^' "*) 
while the relative width of the peaks reduces. 

A similar behaviour, but at much higher gravity (already the first excited state leaves the "global chaos" region at 
7 — 0.625 and enters the coupled oscillators regime at 7 = 27) and moreover complicated by avoided crossings which 
cause deviations from this pattern at some 7 values, is observed for the other states with low quantum numbers. 

V. CONCLUSIONS, PROSPECTIVES AND AKNOWLEDGEMENTS 

I have conducted the first extensive study of the dynamics of the quantum double pendulum: even if not exahustive 
it has allowed us to observe a close correspondence between classical and quantum structures in phase space in all 
three the classical regimes: from the "free rotors" one at low classical gravity parameter 7 to the "coupled harmonic 
oscillators" regime at high 7, all through the "global chaos" regime for intermediate values of 7. In particular, 
notwithstanding the persistence of some regular Husimi functions in the "global chaos" regime, the Husimi functions 
of most of the states in that region are quite complicated, suggesting that the time evolution of quantum packets might 



simulate rather well the chaotic classical evolution, spreading rapidly over most of the phase space and remaining for 
fairly long times in such a state before eventually collapsing again in a localized packet, as expected from the well 
known analysis of Ref. [1] . The study of the time evolution of suitably placed minimum uncertainty packets will be 
the subject of a forthcoming paper. 

Projections of the Husimi functions on other phase space sections -(Ai, (pi) planes for (p2 ^ and (A2, ip2) planes 
for different values of (pi- and investigation of other combinations of the length and mass parameters / and /i might 
give some interesting insight too, especially when compared to the case studied here. 

I wish to thank G. Mantica and S. Locklin for useful comments and suggestions. 
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FIG. 1. An ideal double pendulum. 
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FIG. 2. Comparison of theoretical (dash) and numerical (full line) energy Vs. level number curves. 
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FIG. 3. Energy surfaces in scaled action variable representation. The plots show I2 versus Ii at constant energy h = 1. The 
radii indicate the areas to be calculated to evaluate the density of states. 



12 




gravity parameter 



FIG. 4. Energy curves. 



13 



a 






■I' .' 


a 






' .' 


* 






.■ 


,-■ P"" 




100- 




." ■" 




«■ 


,■' 




,' ■■ ■ 




■' 


(-'. 




IS. 




- -■ ■ 





4.4 46 43 3.0 i.i i.4 



lU Ii7 L^ U7 u ^17 



lU Ii7 L^ U? U 'L17 



FIG. 5. (Color online) An avoided crossing of two levels, one of which is a principal resonance one. The points at which the 
Husimi functions are calculated are marked as bigger dots on the energy curves shown on the right of the figure. The gravity 
parameter increases from top to bottom. For comparison, the classical Poincare section for 7 — 0.04888 is superimposed on 
the corresponding Husimi function. 
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FIG. 6. (Color online) An avoided crossing of three levels, one of which is a principal resonance one. The points at which the 
Husimi functions are calculated are marked as bigger dots on the energy curves shown on the right of the figure. The gravity 
parameter increases from top to bottom. For comparison, the classical Poincare section for 7 = 0.0597877 is superimposed on 
the corresponding Husimi function. 
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FIG. 7. (Color online) An avoided crossing of two levels, one of which is a principal resonance one. The points at which 
the Husimi functions are calculated are marked as bigger dots on the energy curves shown on top of the figure. The gravity 
parameter increases from top to bottom. For comparison, the classical Poincare section for 7 — 0.04888 is superimposed on 
the corresponding Husimi function. 
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FIG. 8. (Color online) Examples of Husimi functions in the "global chaos" region, a): n = 33, 7 — 5.4, 7 — 0.22445. b): 
n = 36, 7 = 5.4, 7 = 0.21950. c;: n = 127, 7 = 9.8, 7 = 0.15692. dj: n = 118, 7 = 8.1, 7 = 0.14632. e^: n = 117, 7 = 8.1, 
7 = 0.14815. /;: n = 74, 7 = 8.6, 7 = 0.19053. 
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